碳离子在不同材料叶片的多叶光栅上产生的次级粒子研究

在碳离子放射治疗中，碳离子束与治疗头设备和患者身体相互作用产生的次级粒子可以到达患者体内的许多区域，在产生的次级粒子中以中子和$\gamma $射线的产额为最大。在不影响束流配送功能的情况下，减少碳离子放疗中产生的次级中子和$\gamma $射线对于降低放疗后出现的正常组织并发症及二次肿瘤风险有着非常重要的意义。本文利用蒙特卡罗(Monte Carlo)方法模拟计算了被动式束流配送系统下，400 MeV/u碳离子束照射到由不同材料叶片构成的多叶光栅(MLC)形成典型的10 cm×10 cm方形射野时，在水模体中产生的次级中子和$\gamma $射线所沉积的剂量及空间分布等。模拟结果显示:碳离子束通过MLC形成射野后在水模体中产生的次级中子主要分布在水模体的入射端，次级$\gamma $射线较为均匀的分布在整个水模体内，且较多分布在具有展宽Bragg峰(SOBP)射野在水模体中贯穿时的坪区。对于MLC叶片材料的选择，则需根据实际情况对叶片厚度以及次级粒子当量剂量的要求来确定。本文通过模拟研究不同MLC叶片材料产生次级粒子的情况，为被动式束流配送系统中MLC叶片及其他元件的材料选择提供了科学依据。     

超临界压力CO2在水平圆管内流动传热数值分析

采用SST k-w低雷诺数湍流模型对加热条件下超临界压力CO2在内径di=22.14 mm,加热长度Lh=2440 mm水平圆管内三维稳态流动与传热特性进行了数值计算.通过超临界CO2在水平圆管内的流动传热实验数据验证了数值模型的可靠性和准确性.首先,研究了超临界压力CO2在水平圆管内的流动传热特点,基于超临界CO2在类临界温度Tpc处发生类液-类气“相变”的假设,揭示了水平圆管顶母线和底母线区域不同的流动传热行为.然后,分析了热流密度qw和质量流速G对水平圆管内超临界压力CO2流动换热的影响,通过获取流体域内的物性分布、速度分布和湍流分布等详细信息,重点解释了不同热流密度qw和质量流速G下顶母线内壁温度Tw,i分布产生差异的传热机理,分析结果确定了类气膜厚度d、类气膜性质、轴向速度u和湍动能k是影响顶母线壁温分布差异的主要因素.研究结果可以为超临界压力CO2换热装置的优化设计和安全运行提供理论指导.     

MoS2/SnS异质结太阳能电池的模拟研究

硫化亚锡(SnS)是一种Ⅳ-Ⅵ族层状化合物半导体材料,其禁带宽度与太阳能电池最佳带隙1.5 eV非常接近,并且在可见光范围内光的吸收系数很大(α>10⁴ cm^-1),因此SnS是一种很有应用前景的材料。本文利用太阳能电池模拟软件wxAMPS模拟了MoS₂/SnS异质结太阳能电池,主要研究SnS吸收层的厚度、掺杂浓度和缺陷态等因素对太阳能电池性能的影响。研究发现:SnS吸收层最佳厚度为2 μm,最佳掺杂浓度为1.0×10¹⁵ cm^-3;同时高斯缺陷态浓度超过1.0×10¹⁵ cm^-3时,电池各项性能参数随着浓度的增加而减小,而带尾缺陷态超过1.0×10¹⁹ cm^-3·eV^-1时,电池性能才开始下降;其中界面缺陷态对太阳能电池影响比较严重,界面缺陷态浓度超过1.0×10¹² cm^-2时,开路电压、短路电流、填充因子和转换效率迅速下降。另外,通过模拟获得的转换效率高达24.87%,开路电压为0.88 V,短路电流为33.4 mA/cm²。由此可知,MoS₂/SnS异质结太阳能电池是一种很有发展潜力的光伏器件结构。     

基于NPP-VIIRS夜间灯光数据的山西省GDP空间化模拟

选取山西省为研究对象，以美国国家极轨合作仪件-可见红外成像辐射计套件（NPP-VIIRS）夜间灯光数据、GDP统计数据等为数据源，构建GDP空间化拟合模型，建立山西省GDP密度图，据此研究山西省经济的空间差异性。通过对NPP-VIIRS夜间灯光数据的空间化处理，提取灯光指数，并将其与GDP进行回归拟合，建立最佳回归模型，得到GDP密度拟合图；利用县级GDP数据进行线性纠正，从而提高GDP的模拟精度。结果表明：（1）NPP-VIIRS夜间灯光数据与GDP的相关性较高，可用于山西省GDP模拟；（2）与GDP分区建模相比，GDP整体建模的精度更高；（3）山西省GDP的空间分布整体呈由城市中心逐渐向周边辐射的特点，构成GDP过渡带。     

一类分数阶Langevin方程block-by-block算法的数值分析

分数阶Langevin方程有重要的科学意义和工程应用价值，基于经典block-by-block算法，求解了一类含有Caputo导数的分数阶Langevin方程的数值解.Block-by-block算法通过引入二次Lagrange基函数插值，构造出逐块收敛的非线性方程组，通过在每一块耦合求得分数阶Langevin方程的数值解.在0<α<1条件下，应用随机Taylor展开证明block-by-block算法是3+α阶收敛的，数值试验表明在不同α和时间步长h取值下，block-by-block算法具有稳定性和收敛性，克服了现有方法求解分数阶Langevin方程速度慢精度低的缺点，表明block-by-block算法求解分数阶Langevin方程是高效的.     

水化硅酸钙纳米压痕分子动力学模拟方法优化

针对水化硅酸钙纳米压痕模型忽视了压头与基底之间相互作用的问题，由尺寸差异引起的金刚石压头难以计算的问题，以及Wittmann模型无法得到实际接触面积的问题，提出了新的模型与计算方法．结合分子动力学方法，采用金刚石压头-Wittmann模型基底的组合方式构建无定形态水化硅酸钙纳米压痕试验模型．在建模阶段，考虑到压头模型与基底模型粒子间尺寸差异，提出了等比例替换模型，通过公式推导并就不同尺寸模拟结果验证了等比例替换模型的可行性．在计算阶段，提出了局部前处理的弛豫方法进行模拟．确定最大荷载位置处的接触面积为546 nm2，进而求出水化硅酸钙模型硬度H为0.84 GPa、折合模量Er为30.52 GPa．并通过纳米压痕试验，验证了模拟结果的准确性，证明了模型的科学性，对今后水化硅酸钙(C-S-H)纳米层面的模拟具有重要借鉴意义．     

峨山褐煤的分子结构和分子模拟

作为重要的化石能源，褐煤资源潜力巨大、分布广泛但综合利用率低。研究褐煤的分子结构模型，有助于预测褐煤在热解、液化和气化过程中的化学反应机理及反应路径，进而提高褐煤的综合应用水平。以云南峨山褐煤为研究对象，利用傅里叶变换红外光谱、13C核磁共振波谱及X射线光电子能谱等分析测试方法，获取了峨山褐煤的含碳、含氧及含氮结构参数。在此基础上，借助Gaussian 09计算平台，采用量子化学建模的方法构建并优化了峨山褐煤的分子结构模型。研究结果表明：峨山褐煤的芳碳率为39.20%，芳香碳结构主要为苯和萘，且芳香桥头碳与周边碳的比值χ_b为0.07；脂碳率为49.51%，脂肪碳结构主要为亚甲基，季碳和氧接脂碳；氧原子主要存在于羟基、醚氧、羰基和羧基结构中；含氮结构则以吡啶为主。基于元素分析、13C 核磁共振波谱分析，又经过热重实验消除褐煤中残余水分的影响后，计算出峨山褐煤的分子式为C₁₅₃H₁₃₇O₃₅N₂。依据分子式及分析结果计算出峨山褐煤的结构单元含量并构建出其初始结构模型，采用半经验法PM 3基组及密度泛函理论M06-2X/3-21G基组对初始分子构型进行优化。优化后的分子模型具有明显的三维立体特征，芳香环之间较为分散且在空间中排列不规则，芳香簇主要通过亚甲基、醚氧基、羰基、酯基和脂肪环连接，含氧官能团主要分布在分子边缘，脂肪族侧链较多。对优化后的分子模型进行振动频率计算进而获得了分子模型的模拟红外光谱，其与实验红外谱图吻合度良好，证明了峨山褐煤分子结构模型的准确性、合理性。分子结构模型的构建有利于直观地了解峨山褐煤的分子结构特征，从而有助于从微观分子角度研究峨山褐煤的宏观性质。同时，峨山褐煤分子结构模型可为其在热解、液化和气化等领域研究中提供理论指导。     

染料敏化太阳能电池载流子传输的数值模拟

由于在染料敏化太阳能电池(dye-sensitized solar cell, DSSC)中存在染料弛豫、半导体薄膜中电子与氧化态染料分子发生反应和电子在电解质中与氧化态离子复合等不利反应,利用一个更完善的DSSC载流子传输模型对电池的光电性能进行模拟就显得非常重要。为此,本文基于由多重俘获理论建立的DSSC中的包括电子、染料阳离子、碘化物和三碘化物在内的载流子传输模型,数值模拟得到了不同TiO₂薄膜厚度、不同入射光强度与不同染料分子吸收系数下DSSC的J-V曲线。结果表明,随着TiO₂薄膜厚度的增加,太阳能电池的短路电流密度增大,开路电压减小,光电转换效率先增大后减小。当DSSC的TiO₂薄膜厚度为20 μm时,光电转换效率达到最大值7.41%,同时光电转换效率随入射光强度与染料分子吸收系数的增大均有一定程度提高,其中在吸收系数为4 500 cm^-1时,光电转换效率为6.73%。以上结果可以为改进DSSC的光电性能提供理论指导。     

A filter method for nonlinear semidefinite programming with global convergence

In this study, a new filter algorithm is presented for solving the nonlinear semidefinite programming. This algorithm is inspired by the classical sequential quadratic programming method. Unlike the traditional filter methods, the sufficient descent is ensured by changing the step size instead of the trust region radius. Under some suitable conditions, the global convergence is obtained. In the end, some numerical experiments are given to show that the algorithm is effective.